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O I Abstract 

' The frame and scale dependence of the pair-term contribution to the electromag- 

netic form factor of a spin-zero composite system of two-fermions is studied within the 
O \ Light Front. The form factor is evaluated from the plus-component of the current in 

the Breit frame, using for the first time a nonconstant, symmetric ansatz for the Bethe- 
Salpeter amplitude. The frame dependence is analyzed by allowing a nonvanishing plus 
component of the momentum transfer, while the dynamical scale is set by the masses of 
^ ' the constituents and by mass and size of the composite system. A transverse momen- 

■ tum distribution, associated with the Bethe-Salpeter amplitude, is introduced which 

allows to define strongly and weakly relativistic systems. In particular, for strongly rel- 
ativistic systems, the pair term vanishes for the Drell-Yan condition, while is dominant 
for momentum transfer along the light-front direction. For a weakly relativistic system, 
fitted to the deuteron scale, the pair term is negligible up to momentum transfers of 
l(GeV/c)^. A comparison with results obtained within the Front-Form Hamiltonian 
dynamics with a fixed number of constituents is also presented. 



To appear in Nucl. Phys. A 



1 Introduction 



Present knowledge of the structure of hadrons and nuclei mainly comes from electroweak 
form factors, in elastic and transition regimes, and from deep-inelastic structure functions. 
To perform a meaningful comparison between theoretical models and the experimental data, 
one needs a description of the bound system of interacting constituents and a consistent 
current operator. Within a Hamiltonian approach, the state of the system is defined on 
a specified hypersurface of the space-time that does not contain timelike directions. Dirac 
identified three space-time hypersurfaces, adequate to define the state of a relativistic system, 
which correspond to different forms of relativistic Hamiltonian Dynamics, namely the Instant 
Form, the Front Form and the Point Form [||]. 

In the Front-Form dynamics, the consistency between the current operator and the 
state of the hadron system has been discussed both from a field theoretical point of view 
and within approaches with a fixed number of particles. In field theory, the state has an 
infinite number of components in the Fock space0. However, for practical applications 
only the lowest Fock component, or valence component, is usually modeled and used in the 
calculations of electroweak form factors. In principle, the infinite set of coupled eigenvalue 
equations for the Hamiltonian operator in the Fock space can be replaced by an effective 
squared mass operator acting in the valence sector; at the same time, it is possible to express 
systematically the higher Fock-state components of the wave function as functionals of the 
lower ones 0, 0, The effective electroweak current operator to be used with the valence 
component of the state can be consistently derived within the framework of the Bethe- 
Salpeter equation projected at equal light-front time, as recently shown in Refs.[|], ||. (For 
recentinvestigation on other aspects of the Bethe-Salpeter equation within the light front, 
see, e.g. 0,1). 

If a fixed number of interacting particles is assumed, then a Front-Form Hamiltonian 
dynamics (FFHD) can be developed. One can use the Bakamjian- Thomas construction!^, 
where an effective interaction in the mass operator can be chosen such that a unitary rep- 
resentation of the Poincare group is possible. Taking the Front-Form spin operator as the 
free one and an interaction in the mass operator that commutes with such a spin operator, 
one is able to explicitly construct generators of the Poincare group||10|, |lT|. In this case, the 
eigenfunctions of the mass operator are normalized to one, differently from the field theoret- 
ical valence component, which has in general a probability less than one[0. A wide number 
of papers have been devoted to the evaluation of properties of hadrons and nuclei within 
the Front-Form Hamiltonian djTiamics with a finite number of particles (see, e.g. Refs. 



T2| , [T3|, |T^, |T^, |T^, just to give a short account of some previous works). In partic- 



ular, it should be noted that FFHD for the two-nucleon case yields the possibility to retain 
the large amount of successful phenomenology developed within a nonrelativistic approxi- 
mation (see, e.g., |T2|, |I^]). Usually, to obtain electroweak form factors of nuclear systems or 



hadrons, the Drell-Yan condition on the momentum transfer, i.e. = + = 0, and the 
matrix elements of the plus component of the current, j"*", are adopted. The physical argu- 
ment often advocated in favour of these assumptions is that the production of pairs from the 
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incoming photon (nonvalence or pair-term contribution) is suppressed in the Drell-Yan frame 
by hght-front momentum conservation (see, e.g., [0). This was indeed proved in schematic 
covariant field theoretical models, for spin-zero two-boson composite systems|2l| and for the 
pion with pseudo-scalar coupling to the quarks||2^. However, for spin-one systems the pair 
term survives and contributes to in the Drell-Yan frame [p3| , p^ , |25[| . Furthermore, in 
this case the pair term is necessary to keep the rotational invariance of the form factors for 
spin-one particles. Within FFHD, in the Drell-Yan frame, the matrix elements of the plus 
component of the current should fulfill the so called angular condition [^]. This constraint is 
not satisfied by a calculation considering only the Front-Form wave function and the impulse 
approximation. This approximation implies an ambiguity in extracting form factors for spin- 
one systems, as shown for the deuteronp6|, 0, ^ and for the rho-meson I^H], | 



30 . To 



solve this problem, related to the lacking of covariance for the impulse approximation current 
operator, some physically inspired combinations of matrix elements were proposed to extract 
form factors from the current, as the "good-component" approach]^ or the elimination of 



spurious contributions to the form factors [p2 



To avoid the difficulties associated with these ambiguities, another approach was 

In these works, it was shown 
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proposed in Ref. ||33| and applied to the deuteron|T9|, 
that a current operator which satisfies the requirements of Poincare, parity and time-reversal 
covariance, as well as hermiticity and current conservation, can be obtained from a one-body 
operator in the Breit frame with momentum transfer along the z-direction. Then, for any 
hadron system the electromagnetic form factors can be calculated from the wave function 
without ambiguities. 

This major development in the calculation of the form factors from the Front-Form 
wave function is, however, confronted with the fact that, even in simple field-theoretical 
models, the pair production mechanism contributes in the ^ Breit frame, as was 
discussed in the calculation of the pion form factor from j"*" in Ref. [^. Therefore, within 
an approach with a fixed number of particles such a mechanism should be taken into account 
through an effective two-body current. 

The aim of the present work is to investigate in a covariant model, based on a non- 
constant, symmetric vertex function, the effect of the pair term in the evaluation of elec- 
tromagnetic form factors for a pion-like system composed by two identical fermions. It is 
well known that, to avoid divergences in the evaluation of the covariant triangle-diagram for 

a pion-gg 



the form factor some regularization has to be introduced. In Refs. 
vertex function, non symmetric in the four-momenta of the quarks, was adopted. However, 
a non symmetric vertex cannot be considered a realistic approximation of a qq bound-state 
amplitude and phenomenological problems arise (e. g., the form factor and the weak decay 
constant cannot be simultaneously reproduced) |^. Another approach has been proposed 
in Ref. 
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where the fermion loop was regulated by considering a non-local photon ver- 
tex. In the present work, for describing the momentum part of the coupling between the 
constituents and the spin-0 system, we use a covariant model with the following form of the 
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vertex function 

C C 
" (A;2 - m| + ^e) ^ ((P - A;)2 - m| + te) ' 



Differently from |22, 23, 24], we adopt a vertex function which is symmetric by the exchange 
of the momentum of the two fermions and imphes a hght-front valence wave function with 
the same property, as shown in the following Sections. In this way, we are simulating the 
symmetry properties of a Bethe-Salpeter amplitude derived from quantum field theory. The 
other main ingredient of our covariant model, for calculating the form factor of the composite 
system, is the electromagnetic current, that is taken in impulse approximation. It should be 
pointed out that within our approach such a current is conserved (see Sect. II). In the spirit 
of Ref. |3^, we study the importance of the pair diagram evaluated in different Breit frames, 
which differ for the direction of the spatial part of the momentum transfer. We analyze and 
compare two systems that have relativistic or nonrelativistic nature, respectively. In order to 
better define this feature, we first construct the valence wave function from the symmetric 
ansatz for the Bethe-Salpeter amplitude and then, from this wave function, we build up 
the transverse-momentum distribution of the constituents. This momentum distribution 
plays a twofold role, allowing one: i) to make contact betwen our covariant model and 
dynamical models of the composite system, developed within approaches with a fixed number 
of particles, and, more important, ii) to quantitatively define the two limiting cases that we 
will consider, namely the strongly and weakly relativistic systems. In the first case, we make 
calculations for a spin-0 model, well suited for the pion. Then we compare these results 
with the ones obtained with a realistic pion wave function generated by a potential able to 
describe the meson spectroscopy]^. In the second case, we artificially adjust the transverse 
momentum distribution of the model to the deuteron scale, to get insight into the pair-term 
contribution for a weakly relativistic system with a small average transverse momentum for 
the constituents (getting rid of the lengthy algebra of the mi ^ 1112 case, corresponding to 
the actual case of charged heavy mesons) . 

The paper is organized as follows. In Sec. II, our model for the spin-zero, two- 
fermion system with a symmetric vertex function is presented. In Sec. Ill, the light-front 
valence wave function for the covariant model is introduced, as well as the corresponding 
transverse momentum distribution and the corresponding elastic form factor. In Sec. IV, 
numerical results are presented for i) the electromagnetic form factor of our model, ii) the 
separate contribution of the pair term and iii) the form factor corresponding to the valence 
component. In Sec. V, we draw our conclusions. 



2 Electromagnetic form factor of a pion-like system 

In our model, the electromagnetic current for a two-fermion composite system with spin 
equal to - i.e., a pion-like system considered as qq bound state - is calculated in one- 
loop approximation (triangle diagram), modelling the Bethe-Salpeter amplitude through a 
symmetric vertex function in momentum space with a pseudoscalar coupling between pion 
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and quark degrees of freedom. This coupling is suggested by a simple effective Lagrangian 
(see, e.g. 0) 

Ci = -ig^ ■ q-f^fq . (2) 

The coupling constant g is given by the Goldberg- Treiman relation at the quark level, namely 
g = m/f.,^, with m the mass of the constituents and /^r the pion decay constant. 

The electromagnetic current of tt"*" is obtained from the covariant expression corre- 
sponding to the triangle diagram (see, e.g., |^ and |^): 



f = -^2e^iV, J ^Jr[S{kh'S{k - P')rS{k - Ph']A{k, P')A{k, P) , (3) 

where Sip) = , A^^c = 3 is the number of colors, P^ and P'^ = P^ + q^ are the 

p — m + le 

initial and final momenta of the system, g'^ is the momentum transfer and k^^ the spectator 
quark momentum. The factor 2 stems from isospin algebra. (Current conservation can be 
easily proven in the Breit frame: after performing the trace in q ■ j, one notices that the 
integrand of the resulting expression is odd by changing k —>■ —k; this means that q ■ j is 
zero.) 

In Eq. (^, we introduce the symmetric vertex function of Eq. (0). This vertex 
function produces a light-front wave function symmetric by the interchange of quark and 
antiquark momenta, and is not affected by the conceptual difficulties associated with the use 
of the nonsymmetric regulator, as mentioned in the Introduction (see Refs. ||30|, p6|). Since 
we have not specified the dynamics which drives the Bethe-Salpeter amplitude, we have to 
resort to a physical condition for normalizing the Bethe-Salpeter vertex. As a matter of fact, 
the normalization constant C in the vertex function, Eq. (|I]), is fixed by imposing the charge 
normalization condition (i.e. the pion form factor at zero momentum transfer must be equal 
to 1). 

In our analysis we consider Breit frames, where the momentum transfer has the 
spatial component parallel to the z — x plane. By using Front-Form variables, i.e. k~^ = 
k^ + k^ , k~ = k^ — k^ , k± = {k^, k"^) one has 



q = —q = y— g^sina, q^ = y— g^cosa, g^ = 
The Drell-Yan condition g^ = is recovered with a = 0, while the q'^ = \J ~q^ condition [^| 



comes with a = 90°. (Note that the angle 6 of Ref.p6| corresponds to a + 90°). The initial 



and final momenta of the composite spin-0 bound state are: P^ = E = E' = — g^/4, 

n = -P± = gl/2 and Pi = -P, = g+/2. 

As well known, the pion form factor can be extracted from the covariant expression: 

j'^ = e(P^ + P'^)P^(g2) . (5) 
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If covariance and current conservation are fulfilled in a model calculation, one can obviously 
use any frame and any nonvanisliing component of the current to calculate the electromag- 
netic form factor. We calculate the form factor for our pion model defined by Eq. (|l|) and 
Eq. using the plus component of the current in the Breit frame with q in the z — x plane. 
In the evaluation of the form factor one can single out two nonvanishing contributions in Eq. 
i) iH, 1111,11: 

FM') = Fi'\q\a) + FP{q\a) , (6) 

where Fy\q'^,a) has the loop integration on constrained by < < P^, see the 
light-front time-ordered diagram in Fig. 1(a), while F^^^\q'^,a) has the loop integration 
on k~^ in the interval < k^ < P"*", see Fig. 1(b). The valence component of the 
pion contributes to Fy\q'^,a) only, but in our model it does not give the full result in 
< k~^ < P~^, as discussed in detail in Sec. IV. The component Fy^\q'^,a) of the form 
factor is the contribution of the pair production mechanism from an incoming virtual photon 
with g+ > 0. 

The two contributions to the form factor obtained from j"*" are given by the following 
expressions 



2 



N, f (Pk^dk+dk-e{k+)e{P+ - A;^ 



(P+ + P'+)/2 (27r)4 J k+{P+ - k+){P'+ - A;^ 
and 

_ rr? r d'k^dk^dk-Qjk^ - P^P'^ - k^) 

^^'"^ - '\p^ + p'^)PA2^y] fc+(p+-fc+)(p'+-fc+) n(/c,p,p(^) 

where 

Tr[0+]k{k,P)K{k,P') 



n(fc,p,p') 



(fc- - fc- + te){P- -k--iP- /c)- + 

1 



{P'- -k^- {P' - k)-n + *e) ' 
with the on-energy-shell values of the individual momenta given by 



(9) 



^ , (P - ^ ^^^^ and (P' - ^ . (10, 

In Eq. (y), the trace Tr[(9+] of the operator 

0+ = {f + m)7^(|^ -P' + m)-f+{^ -p + m)j^ , (11) 

is given by: 

lrr[C+] = -k-{P'+ -k+){P+ -k+) + {kl + m^){k^ - P+ - P'+) 

- \k^-{P'^-P^){P'^-P^) + \kUl. (12) 
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The explicit form of the symmetric regulator function in Front-Form momentum co- 
ordinates is given by 



A{k,P) = C 



+ C 



k-+{k- 



k'j_ + m\ — It ^ 



k+ 



k+){p- -k- - 



k)l + ml 



(13) 



where the position of the poles for k~ clearly appears. 

The detailed forms of F'^^ and Fy^\ after integrating over are given in the 
Appendices A and B, respectively. In what follows, we will discuss some general features of 
Eq. (0) and Eq. (|). 

Since the integration range of k^ is < k^ < P^ in Eq. (|^) and P+ < k^ < P ^ 
in Eq. then, the sign of the imaginary part of some of the poles in the ^"-complex plane 
changes (see Eq. (^ and Eq. (p^). The poles that have their imaginary part modified are 



^(1) — P {P ^)on + 



P+ - k+ 



(P - k)l + 
P+ - k+ 



P+ - A;^ 



and 



^(2) — P 



(P - k)\ + m 



R 



+ 



le 



P+-k+ P+-k+ 



(14) 



(15) 



The last one comes from the vertex function, Eq. ([131) . The difference in the sign of the 
imaginary parts of fc^^ and A;^^ for the intervals < A;^ < P"*" and P"*" < fc"*" < P 
is the mathematical signature of the pair production mechanism, which appears just in the 
second interval. 

The sum of the contributions Py^(g^,a) and Fy^\q'^,a) yields the covariant result, 
dependent upon only. Then the different directions of q in the Breit frame can only change 
the values of F^^q"^, a) and Fj^^^\q'^, a), but not their sum. For instance, by chosing g"*" = 
(i.e. a = 0) Fj^^^\q'^,a) vanishes and therefore F^^\q'^,a) alone gives the whole, covariant 
result 

It is interesting to note that in our model the pair term is linear in g+, for small q^, 
as one can verify by direct inspection of the structure of Eq. (|]), once the k~ integration is 
performed. As a matter of fact, the contour for the Cauchy integration in the calculation of 
the pair diagram can be closed in the upper complex k~ semi-plane (see Appendix B), and 
consequently the poles in the integrand of Eq. (|) are 



^(3) 



P'- - (P' - k)-^ + ze = P' 



(P' 



A;)5_ + 



P'+ - k+ 



and 



-(4) 



P'- 



(P' 



P'+ - k+ 



(16) 



(17) 
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The first pole, Eq. (0), comes from tlie last factor of Eq. (||), while the second pole, Eq. 
([T7| ) comes from A{P',k). The position of both poles in k~ are ~ {q~^)~^ in the limit of 
— >■ 0. Then, in order to find the dependence of F^^^\q^,a) on q'^ in this limit, it is 
enough to count the power of g"*" in Eq. (H), when the residues are evaluated. The phase- 
space factor in the denominator of Eq. (H) is of the order of {q^Y- Then, let us consider the 
trace, Eq. (|12|) , with the proper values of k~ (see Appendix B). The first term and third one 
are ~ g"*", the second term and the fourth one are of order {q'^Y- Then the trace is of order 
{q^Y ■ Evaluating the contribution to the residues from the remaining part of Eq. (P), we 
found that n(/c, P, P') is of the order of (g"*")^. Therefore the integrand goes to a constant 
for g+ 0, and thus F^^^\q^^ a) is proportional to q'^ in this limit, because of the range of 
the integration. 

A relevant feature in the analysis of the form factor is given by the presence of a 
contribution which is instantaneous in the light-front time, and is produced by the instanta- 
neous term present in the Dirac propagator. As a matter of fact, the Dirac propagator can 
be decomposed using the Front-Form momenta as follows 



k^-m^ + te k+{k- -k-^ + ^) 2k 

where the second term, proportional to 7^, is an instantaneous term in the light-front time. 
It should be pointed out that the instantaneous contribution to the form factor is produced 
only by the spectator fermion. Indeed, the instantaneous terms pertaining to the other 
propagators do not contribute, because of the factor 7+ from the current and the property 
(7+)^ = 0. 



In our symmetric model, the instantaneous term of Eq. (p!8D contributes both to 
Fy\q'^,a) and F^^^\q'^,a), due to the analytic structure of the symmetric vertex function 
of Eq. (|l]). These contributions are of nonvalence nature, since they cannot be reduced to 
the impulse approximation with the valence wave function. 



3 Valence Light- front wave function 

The valence component of the light-front wave function can be obtained from the Bethe- 
Salpeter amplitude eliminating the relative light-front time, i.e., constraining to equal light- 
front time the external space-time coordinates of the two fermions, after dropping the in- 
stantaneous terms of the external Dirac propagators^. Physically, in the external legs only 
light-front propagating particles are allowed. It is worth noting that the effect of the instan- 
taneous terms, which are present in a Bethe-Salpeter approach, could be included in effective 
(many-body) operators to be used together with the valence wave function ||^. 
In the present model the Bethe-Salpeter amplitude is 

m k + m c:.,, , k — P + m , , 

^ fc, P = ^ ^\ 7 A A;, P „ . 19 
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— * — * 

The momentum part of the valence component of the hght-front wave function, k±; P+, Pi 

can be obtained by ehminating out from Eq. (|T^: i) the instantaneous terms, ii) the factors 
containing gamma matrices in the numerator, and iii) the and (P"*" — fc"*") factors appear- 
ing in the denominator. Then, after introducing the exphcit expression for A, one has to 
integrate over k~ , viz. 



dk- 



27T {k- 

1 



^on+fT)(-P ^ (P ^)on + p+-k+ . 



X 



+ 



where A/" is a normahzation factor 



(P — ky — m|j + 



(20) 



(21) 



Performing the k integration in Eq. (|20|) , one has 

$(A;+,^^;P+,P^) 



Mi 



Af 



< 1; M^imlml) 



{P-k)\+ml 



+ 



1-x 



{22) 



Pf ; and 



where x = k~^/P^, with < x 

the square of the free mass is Mq = A^^(m^,m^). Since we have chosen a nonconstant, 
symmetric A, a second term appears in Eq. 



P^) , differently from Ref. p2|, and then the 
momentum part of the wave function becomes symmetric by the exchange of the momenta 
of the two constituents. 

By using only the valence component and generalizing the results of Refs. 



22 



the electromagnetic form factor, P^^^\ evaluated in the Breit frame is written as follows 



a) 



2n^{P'+ + PH 



X 



X 



(fk^dk+e{k+)9{p+ - k+) 
k+{p+ - k+){p'+ - k+) 



$(fc+,fcx;P'+,^) 



/+ 



^{k+,kr.P^ 



2 



--) 

2 ' ' 



P^ 



4 



(23) 



with k~^ = {k'j_ + irP')/k~^ (see Eq. (p!0|)). Once the normalization constant C is determined 
from the condition P7r(0) = 1 in Eq. (|^), the value of P^'^-^) for = yields the probability 
(independent of a) of the valence qq component in the pion, rj. It should be pointed out 
that a value of ?7 < 1 is expected, if the non- valence contributions are important, and this 
is just what occurs in our model (see Sect. IV for details), differently from the case of a 
nonsymmetric ansatz for the Bethe-Salpeter amplitude, where r/ = 1 p2| . 
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Equation (p3D represents a point of contact between a field theoretical approach and 
the FFHD with a fixed number of particles, adopting an impulse approximation current 
operator. Indeed, including a proper factor in $, namely Mq / P"*" and normalizing the wave 
function one can recover the FFHD expression for the form factor, p^^^^ _ in particular, 
putting a = 0°, one obtains the standard FFHD expression in the frame g"*" = |T3| . 

The valence component of the light-front wave function, Eq. (pOD, is not an eigen- 
function of the total angular momentum, since it is only one of the components of the pion 
state in the Fock space. Therefore, we cannot directly compare the valence wave function of 
the present model with the wave functions corresponding to approaches with a fixed number 

However, 
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of particles IT^, which are eigenfunctions of the intrinsic angular momentum 
in order to make contact with dynamical models, we introduce the transverse momentum 
probability density 



-L 



4vr3 



2lT 



dk+M^ 



k^{mT, — A;+) 



$2(A;+,^^;m,,0) 



(24) 



By integrating f{k±) over k_i_, we obtain the probability of the valence component in the 
pion: 



V 



dk±k±f{kj 



(25) 



Furthermore, the transverse momentum probability density vs fc^/m results to be very 
useful for a quantitative definition of strongly and weakly relativistic composite systems. 
Indeed high values of the transverse momentum distribution for k±/m > 1 are a distinctive 
feature of a strongly relativistic system. In particular, in Sect. IV f{k±) will allow us to 
investigate the influence of the dynamical scale of the system on the role played by the pair 
term. 

Another relevant quantity to be used for constraining the parameters of our model 
is the pion decay constant, /yr. It is defined through the matrix element of the partially 
conserved axial-vector current 



P^<0|Af|7r,- >=imlfJi^ 



(26) 



Following Ref.[^, we take Af = q'y^'y^^q and adopt our ansatz for the pion-gg vertex 
function. In this way we obtain 



m 



I ■i^ATrWS{k)^'S{k - P)]Aik, P) 



and integrating on one gets in terms of the valence component of our model: 



(27) 



4v]-3 



d'^k±dk~^ 
k+im^ — k' 



-$(A;+,fcx;"^^,0) . 



(2J 
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4 Numerical Results 



4.1 Pion model 



In our model, we have two free parameters for the pion: the constituent quark mass, m, 
and the regulator mass, m/j. The constituent quark mass is chosen as m = 0.220 GeV, 
adequate for the meson phenomenology] 37, ^ The regulator mass m^j = 0.6 GeV is 
found by fitting Eq. to the experimental value /^^^ = 92.4 MeV. For the pion mass we 



use the experimental value of 0.140 GeV. As a consequence, the charge radius, obtained from 
(r^) = G^i^TT, comes out to be 0.74 fm, which is about 10% larger than the experimental 



value 



dq 

r' 

exp 



0.67 ±0.02 fm 



In Fig. 2, the results for the pion form factor are shown and compared to the data 



of Refs. p^, The full- model calculations, Eq. (g), nicely agree with the 

new data for the pion form factor |44]. Therefore, our model, based on a non-constant. 



symmetric vertex can reproduce the form factor data consistently with the experimental 
value of /tt, while for the nonsymmetric regulator this was not possible]^]. Remarkably, in 
order to reproduce simultaneously /^r and the experimental form factor within our model, the 
constituent quark mass has to be chosen in the range between 0.2 and 0.3 GeV. Compatibility 
between form factor and decay constant has been also achieved in [RH], where a constant 
pion-gg vertex and a nonlocal photon vertex were adopted. 

The form factor in the Breit frame with a = 0° {q~^ = 0), where the pair-term contri- 
bution is zero, is identical to the sum of F^^ and Fy^\ calculated in the frame corresponding 
to a = 90°, as it must be for a covariant model. In Fig. 2, and Fy^\ for a = 90°, are 
also shown. Differently from the case a = 0°, for a = 90° the form factor is dominated by 
the pair production process, except near = 0. At high values of the momentum transfer 
the form factor is completely exhausted by the pair-term contribution. It is worth noting 
that such a dominance is mainly due to a kinematical effect and appears to be fairly model 
independent (see also Ref.[^). A qualitative argument, which is applicable to the reference 
frame with g"*" = \/ —q^ and essentially follows the work of Sawickipl] 
follows. 

Since the form factor is dimensionless, one can write qualitatively 

dk+ P+ 



is given in what 



10 P+ -I- P'+ 
At the same time, see the discussion in Sect. II, 
amplitude to the form factor, is roughly 

rP'+ dk+ 



P>.^P.- (29) 
the contribution of the pair production 



a) 



(30) 



lp+ P+ + P'+ P'+ + P+ 
The validity of such an approximation is related to the momentum cutoff in the vertex 
function, and therefore it is not reliable for large values of q^. Recalling that for a = 90° 



ml + 



q-^ 



q- 



(31) 
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the estimate of the ratio of the above contributions is given by 



Fy)(g^a = 90°) ^ 




(32) 



At the quahtative level, we can roughly say from Eq. (|3^) that the two contributions to 
the form factor are expected to have about the same magnitude when — = —qfj/u-j = 
ml/2, which gives for the pion —qfj/jj^ = 0.01 (GeV/c)^. Our model calculation yields 
-Qfi/ii) = 0.03 (GeV/c)^. In Ref. ||3^, for the kaon it was shown that the pair term contri- 



bution becomes dominant for —qfj/jj-j ~ 0.2 (GeV/c) , while our estimate (Eq. (|32D) gives 
—qfj/ii) = = 0.13 (GeV/c)^. Therefore, one could argue that our crude estimate is 

able to give a momentum scale where sizeable effects due to the pair production term are 
expected for a = 90°. 

The contributions of the instantaneous part of the Dirac propagator to Fy\q'^,a) 
and Fy^)(g^,a), called F^^l^^^lq"^ , a) and F^^llg^{q'^,a), respectively, are also shown in Fig. 
2 for a = 90°. One can physically understand why in F^^^\q^,a) the instantaneous part is 
important and dominates at high momentum transfers. The interpretation is the following: 
in principle, see Fig. 1 (b), the spectator quark can be exchanged between incoming and 
outgoing pion at a given instant x"*", while the pair of quark and antiquark has been produced 
by the virtual photon at an earlier stage. As the magnitude of the momentum q~{= —q^) 
increases, the time fluctuation for the virtual process decreases and favours the instantaneous 
exchange of the spectator quark between the initial and final pion. In fact. Fig. 2 shows 
the dominance of Fj^^H^^lq'^, a) in the pion form factor as the momentum transfer increases. 
It is worth noting that the value of Fy]„^^(g^,a) is nonzero because of the specific analytic 
structure of the vertex. As a matter of fact, ^^^^^^(g^, a) is nonzero, because of the presence 
of a pole in Eq. (^ at k~ = {k\ + m'j^)/k~^ in the vertex function A, Eq. (|T3D , when one 
choses to close the contour for the Cauchy integration in the lower complex semi-plane of 
k-. 

In Fig. 3, the results for the various contributions to the pion form factor for — g^ = 1 
(GeV/c)^ as a function of the angle a are shown. For increasing angles, the form factor 
changes smoothly from valence to pair-term or nonvalence dominance. 

In Fig. 4, we compare the results for F^), Eq. (g), and F^^^\ Eq. (H), where 
the light-front valence component of the model, defined according to Eq. ( p2l) is used. The 
absolute normalization of the valence component rj, i.e., the probability of the gg Fock-state 
component in the pion, is calculated to he rj = 0.77, differently from the nonsymmetric 



regulator model of Ref. |2^, where r] = 1. The symmetric form of the vertex implies 
contributions from many poles to the form factor and the presence of many poles makes the 
valence component comparatively smaller than in the nonsymmetric case. In a previous work 
on DIS, based on a wave function contribution only, ||3^ a renormalization around 0.5 - 0.75 
was necessary to fit the data. To compare F^^ and Fj^^^\ we have arbitrarily normalized 
p(WF)(^Q^ to 1. As shown in Fig. 4, the momentum behaviour of and F^^^^ is almost 
the same, independently of the reference frame, for a between 0° and 90°. This means that 
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our kinematical argument, about the suppression of Fy\q'^, a) with respect to the full form 
factor in the g+ = \/ —(p' frame, could be extended to F^^^iif' ^ a) as well. 

Once F^^^''(0) is arbitrarily normalized to 1, namely C C / y/r] in Eq. (0), a good 
description of the form factor data for a = 0" is achieved, as one can deduce from Fig. 4, but 
y^H/F) _ ]^gg ]\/[gY^ which overestimates the experimental value, since f2^^^ = fw/ y^- The 
value of f^^^^ is similar to the one found for nonsymmetric models, once the form factor 
below 2 (GeV/c)2 is fitted (see e. g. pi). 



As discussed in Sect. Ill, Fj^^^^ can be formally related to the form factor obtained 
within FFHD, where only the valence contribution appears. Since the momentum behaviour 
of and F^^^) is almost the same, then one can argue that p^^^^ should correspond to 
the contribution of -Fy^ This means that F^^^^ calculated in the Drell-Yan frame could 
represent a good effective approach for evaluating the form factor of the pion. Let us note 
that in FFHD the pair term can appear only as a contribution from two-body currents. 

4.2 Weakly relativistic systems 

In order to investigate the sensitivity of the pair term upon the dynamical scale of the 
composite system, we consider a case sharply different from the pion one. We adjust our 
model to the deuteron scale, and use for the mass of the system the value uid = 1.874 
GeV and for the mass of the constituents m = 0.938 GeV. We adopt a regulator mass of 
mji = 1.1 GeV that gives a mean square radius of 3.25 fm^ comparable to the difference, 
''^D,exp ~ ^p,expy bctwcen the experimental values of the deuteron and proton mean square 
radii. In Fig. 5, we show the results of the form factor calculations for our mock deuteron. 
According to our qualitative kinematical estimate given in the previous subsection (see Eq. 
(p^), with m£, replacing m^), at — ~ 2 (GeV/c)^ the pair term in the Breit frame with 
g+ = \J —q^ is expected to be as large as F^^^ . As a matter of fact, at this momentum 
transfer, the pair term is only about 15% of the form factor, as shown in Fig. 5. Indeed, 
the kinematical estimate does not work for a system with a rapidly decreasing momentum 
distribution. Moreover, we observe that the instantaneous term of the Dirac propagator gives 
a small contribution to the pair term, as well as to F^^\q^^ a = 90°). This fact is related to 
the basically nonrelativistic nature of the costituents in our mock deuteron. As illustrated 
in Fig. 6, the dimensionless product of times the transverse probability density, Eq. (|2^) , 
is peaked at k±_ = 0.06 m for our model. For the sake of comparison, in Fig. 6 it is also 
shown the transverse probability density for an actual deuteron; in particular the density 
has been calculated from the deuteron wave function corresponding to the Reid soft-core 



model of the deuteron |[49||. Remarkably, the same overall behaviour has been obtained in 
both cases. On the contrary the pion is a strongly relativistic system. For our pion model 
k±f{k±) is peaked at k± ~ m, while for the FFHD model of Refs.[0, the peak occurs 



around k± = 2.5m. The FFHD model of Refs.|16, 17] is based on a mass operator with the 



effective interaction of Ref. p7||, which includes one-gluon-exchange and linear confinement 
terms. The pion model from this mass operator, in the g"*" = frame, has a charge radius 
of 0.46 fm for pointlike constituent quarks. This value largely explains why the position of 
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the peak appears at an higher value of the transverse momentum with respect to our model, 
which has a radius of 0.74 fm. 

In summary, for a weakly relativistic system with spin equal to 0, the valence con- 
tribution to the form factor yields a good approximation of the whole form factor up to 
—q^ = 2 (GeV/c)'^ in any reference frame. Therefore, one can argue that the evaluation of 
pFFHD ^^^^ considered a reliable approximation to the full form factor in this case. 



4.3 Pair-term contribution and FFHD estimate 

In order to study in more detail the pair-term contribution, we show in Fig. 7 the ratio 
between F^^^^q"^, a = 90°) and the full form factor. We compare the results of our covariant 
model 

Ri-ov> ^ P""^"^--' «°°) (33) 



with the kinematical estimate, Eqs. (|29|) and (pO|), given by 

= . (34) 

This kinematical estimate of the pair term works fairly well for the pion, while it overesti- 
mates the contribution of F^^^^ for the weakly relativistic system, due to the rapid fall-off of 
the momentum distribution which makes Eq. (|30| ) unreliable. In Fig. 7, we also report the 
following ratio 

where F^^^\q^) and Fj^^^^^q"^) are the pion form factor evaluated in the Drell-Yan frame 
and in the q~^ = \/ —q^ frame, respectively, within the FFHD model of Ref . , using point- 
like constituent quarks and the one-body current. Such a ratio could yield insight on the 
relevance in FFHD of two-body current contributions corresponding to the pair term. In- 
deed, F^^^{(^) represents an estimate of Fy\q'^,a = 0°) = F^{q^), while F^^^^\q^) should 
describe F^^\q^^a = 90°) (cf. comments to Fig. 4) and therefore one can argue that 
^{FFHD) _ ^(coy) _ R{KiN)^ ^ matter of fact, the ratio R^fphd) q^ahtatively agrees 

with the results of our covariant model and with the model-independent kinematical estimate 
(see Fig. 7). Thus, the suppression of F^^^^ in the g"*" = \/—q^ frame with respect to the 
Drell-Yan frame can be explained in terms of a missing contribution in the current, that we 
could identify with two-body current contributions related to the pair term. 

For the sake of completeness, and R^^^^^i are also shown in Fig. 7 for a 

weakly relativistic sytem. Since >> R^'-^'^^\ one can rely only on model calculations 

for inferring the dominance of the valence contribution with respect to the nonvalence, over 
a wide range of momentum transfer. We have also evaluated R^^^^^^ for our mock deuteron 
and we have found R^^^^^^ ~ R^^^^\ This result confirms the small effect of the pair term 
for a spin-0 weakly relativistic system. 
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5 Conclusion 



Within a covariant model with a vertex function which is symmetric in the momentum space 
of the constituents, we have calculated the electromagnetic form factor for a two-fermion, 
pion-like system. Such a model has allowed us to perform a detailed analysis of the pair- 
term contribution both changing the orientation of the Breit frame, where the form factor 
is calculated, and the dynamical scale of the system. 

We have used a pseudoscalar coupling for the pion-gg vertex and a current for the 
point-like constituents in impulse approximation. It is worth noting that in our approach 
the current of the whole system is conserved. We have adopted for the first time a noncon- 
stant, symmetric ansatz for the Bethe-Salpeter amplitude, and we have obtained the valence 
component of the pion state by projecting out the Bethe-Salpeter amplitude on the x"*" = 
hypersurface. The symmetric vertex function implies a valence wave function with the same 
property. In this way we have overcomed some previous conceptual and phenomenologi- 
cal difficulties related to the use of a nonsymmetric regulator for describing the pion. As 
a matter of fact, the form factor and the weak decay constant cannot be simultaneously 



reproduced with a nonsymmetric vertex function (see, e. g., ||22|| ). 

In our covariant model for the pion, the two free parameters, the constituent quark 
mass (m) and the regulator mass {m^), have been fixed as follows: i) the constituent quark 



mass was chosen as m = 0.220 GeV from the meson spectroscopy |37], ii) the regulator mass 



mR = 0.6 GeV was found by fitting the experimental weak decay constant, /^"^^ = 92.4 MeV. 
As a consequence, the form factor obtained in our covariant model nicely agrees with the 
new data for the pion form factor |44]. It should be pointed out that, if we wish to reproduce 
simultaneously and the form factor data with our approach, the constituent quark mass 
should be chosen in the range between 0.2 and 0.3 GeV. Another interesting outcome of our 
symmetric model is that the probability of the pion valence component, r], is less than 1 



(about 75%), at variance with previous covariant calculations p2|, ^ where 77 = 1. The 



em form factor of the pion-like system has been decomposed in two contributions, with the 
second one given by the pair-diagram term. The separation of the covariant result in such 
contributions is unique and does not depend on any particular dynamical model used to 
generate our choice of BS vertex, as long as the four- dimensional impulse approximation is 
adopted to calculate the em current (see Eq. (||) ). 

Following Ref. but with a nonconstant, symmetric vertex function, we have 

investigated the frame dependence of the pair-term contribution to the form factor of pion- 
like systems. Our conclusions strengthen and generalize the conclusions of Ref. In 



particular, after introducing a transverse momentum distribution to distinguish between 
strongly and weakly relativistic systems, we have also analyzed the form factor for the 
second case. 

In general the magnitude of the pair-term contribution in the electromagnetic current 
depends on the frame, because a dynamical transformation change the relative weight of 
valence and nonvalence contributions in the Fock state and then in the form factor, while 
the sum, that is an invariant quantity, does not. For strongly relativistic systems, as the 
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pion, the effect is dramatic: in the Breit frame where = \/ —q^ the form factor is largely 
dominated by the pair diagram for — 7^ 0, while adjusting the model to the deuteron scale 
(weakly relativistic system), the pair term becomes negligible at low momentum transfer and 
only contributes by 30% at -g^ = 4 {GeVjcY. 

The rapid fall-off of the valence component contribution to the form factor of the 
pion for purely longitudinal momentum transfers, = \/—q'^, is a feature which can be 
understood also by a general, crude kinematical estimate. As matter of fact, the ratio 
between the pair contribution and the full form factor is roughly given by (see also, e. g.. 



From such an estimate one can have some insight on the role played by the mass of the 
composite system, in determining the onset of the importance of the pair diagram in the 
form factor. 

Finally, we have compared the results of our covariant model for the pion form factor 
to the ones developed within the Front Form Hamiltonian Dynamics, based on one-body 
current operators. Such a comparison has shown that some general features are shared. In 
particular, the fall-off of the valence contribution in the g+ = \/ —q^ frame faster than the 
one in the g+ = frame can be recovered in the FFHD approach. It should be pointed out 
that the FFHD calculations do not include the contribution of the pair term to the form 
factor, which should appear as a two-body current contribution. 

In the case of the deuteron scale (weakly relativistic system), the relevance of the pair- 
term contribution is shifted towards higher values of the momentum transfer (2-4 (GeV/c)^). 
Therefore, these encouraging results of our exploratory analysis at the deuteron scale urge 
studies which take into account the vector nature of the deuteron. 

The investigation of the pair-term contribution carried out within the covariant model 
of the present work could give a first suggestion for an explicit form of the two-body current 
contribution to be used within phenomenological FFHD approaches. A detailed analysis 
of the light-front pair term in a given field theoretical model can be done using the quasi- 
potential approach of Refs. 0. Within this approach, the light-front pair diagram for 
nonvanishing q^ will come from the contribution of a two-body current, which can be formally 
derived in a consistent way. Numerical investigations of a dynamical model, featuring the 
mentioned properties and suitable for applications to hadrons, will be presented elsewhere. 
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A Analytic Integration on k- in F^) 



In this Appendix, we show in detail how to perform the contour integration on k in Eq. 
(0) for a). This quantity can be rewritten as follows 

Fi'\q\ a) = + F^)^ + Ff)^ + F^)"^ . (36) 



The first term in Eq. (P§), FP", is 



^ ~ 2n J (k+r{P+ - k+){P'+ - k+) ^ ^ 

1 



(k- - k-^ + te)iP- -k--iP- k)-, + le) 

1 



(37) 



(p/- _ fc- - (P' - fc)- + te){k- ~k-^ + tef ' 
where we exphcitly wrote the dependence upon k~ in Tr[0^ {k~)] and 

= , (38) 

The on-energy-shell values of the individual momenta k~^, (P — k)~^ and (P' — k)~^ are given 
by Eq. (p!0|). The normalization factor is 

"^^^^ (39) 
(P+ + P'+)/2(27r)3 ^ ^ 

For the sake of algebraic simplicity the contour of integration in Eq. (^) can be 
closed in the upper complex semi-plane of k~, where only the poles k^j^^ and k'^^^, given by 
Eqs. (|1^) and (|16|), respectively, are present. By evaluating the residues of the integrand in 
Eq. (|37D, one has 



piDa ^ f <fk^dk^d{k+)e{p+ - k- 

J (k+)HP+ -k+)(P'+ -k 



(fc+)3(p+ _ A;+)(p'+ _ k+) {P'- -P- + {P- k)-^ - (P' - k)-J 

Tr[0+{p--{P-k)-J] 



(P- _ (p _ k)-^ _ A;-J(P- - (P - A:)-„ - k-^Y 

Tr[0^{P'- - jP' - fc)-J] 

iP'- - iP' - k)-^ - Kn){P'- - iP' - - kn? 
The second term in Eq. (|36|), Py-**, is given by 

p(/)^ _ z f d^fc^dfc-^rffc-g(fc+)g(P+ - A:+) 

~ 27r y (A;+)2(P+-A;+)(P'+-A;+)2 '"^^^ >^ 

1 



(40) 



(A;- - fc- + ^6)(P- - fc- - (P - fc)-„ + ie){P'~ -k-- (P' - k)-^ + ze) 

(41) 



(A;- -k^ + ie){P'~ - k- - {P' - k)^ + le) 
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where 



P'+ - kr 



(42) 



In Eq. (^TD, the integration contour is closed in the lower complex semi-plane of where 
the poles k~j^^-^ = — le and fc^"^-, = k]^ — le are present. By calculating the residues of the 
integrand in Eq. (^T]), one gets 



-N 



(Pk^dk+e{k+)e{p+ - k+) 
{k+f{p+ - k+)(p'+ - k+y 
1 



X 



X 



(P- _ (p _ A:)-„ - k-JiP'- - iP' - k)-,, - k 
1 

+ 



The third term in Eq. (pG]), F'^^^'^^ is given by 



(43) 



X 



d^k^dk+dk-e{k+)e{p+ - k+) 
27t" J {k+y{p+ - k+){p'+ - k+y 

1 



Tr[0+{k-)] 



(k- - fc- + ie){P- -k--{P- k)-^ + it){P'- -k-- {P' - k)-^ + te) 
1 



(44) 



^k--k'^ + ie){P--k--{P-k)-^ + ie) 
Eq. (^11) is identical to Eq. (^) with P' ^ P. Consequently we can write 

The last term in Eq. (|36|), F^^^'^, that represents the full for the nonsymmetric 
model of Ref . , is given by 



2% 



d'^k^dk+dk'9{k+)e{P+ - A;+ 
k+{P+ -k+f{P'+ -k+f 

1 



-Tr[0+{k-)] 



X 



X 



(fc- - /c- + Z6)(P- - A;- - (P - fc)- + 2e)(P'- - fc- - (P' - /.)" + ze) 

1 

(P- - A;- - (P - A:)^ + ze)(P'- - k' - (P' - A:)^ + le) ' 



where 



(P - k)\ + 
P+-k+ 



(46) 



(47) 
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In Eq. (^61), the integration contour can be closed in the lower complex semi-plane of k~ , 
where only the pole fc^"^-) = — le is present. By evaluating the residue of the integrand in 
Eq. (|46D , one has 

^ - 'V k+iP+ - k+Y{p'+ - k+Y ^ ^ "^'^ 

1 

X 



{P ^on {P ^)on){P' ^on i^' ^)on) 

(P- - Kn -{P- k)n){P'- - k-^ - (P' - k)-^) ■ 
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B Analytic Integration on k- in F^^) 

To perform the contour integration on k~ in Eq. (^) for the pair term, we 
rewrite such a contribution as follows 

F^^) (g^ a) = Fi^'^'^ + Fi^^^' + F^^'^^ + F^^)'^ . (49) 

All the terms in Eq. (^) can be evaluated by closing the integration contour in the upper 
complex semi-plane of k~ . The first term in Eq. (^Of ), F^^^"" is 

" 2n J (k+r(p+ - k+){P'+ - k+) ^ ^ 

1 



(k--k~^ + te)iP--k--iP-ky^ 



1 

iP'- -k-- {P' - k)-^ + ie){k- -k-^ + leY ' 

where the on-energy-shell values of the individual momenta /c~„, (P — /c)~„ and (P' — A;)~„ 
are again given by Eq. (|T0|) and /c]^ is defined by Eq. (^8|) . The normalization factor is the 
same as in Eq. (|39|). 

In Eq. (^Of), only the pole k'^^-^ given by Eq. ( |16|) is present in the upper complex 
semi-plane of k~ . By evaluating the residue of the integrand in Eq. (|49|), one gets 

^ f <fk^dk^d{k+ - p+)e{P'+ - fc+) rr[0+(p-- - (p- - k)-j] 

J (fc+)3(p+_^+)(p'+_A;+) (P'- - (P' - A;)-, - A:-J 

1 

- i^'- -iP- k)-n + - k)-n){P'- - {P' - k)-n ' ■ ^^^^ 
The second term in Eq. F^^^^, is given by 

1 

{k- - k-^ + ie){P- -k--{P- k)-^ - ie){P'- ~k-- (P' - A;)- + le) 
1 

(A;^ - fc^J + 2e)(P'- - A;- - (P' - A;)^ + le) ' ^^^^ 

In Eq. (^2]), the poles A;^-|, Eq. (0), and A;^^-^, Eq. (]T7|), contribute. By calculating the 
residues of the integrand in Eq. ( ^21) , one has 

P(-)^ = -N [ d'^-d^'^^k' - ^")^(^" - ^')Tr[0-iP'- - (P- - k)- )] 
7 (A;+)2(P+- A;+)(P'+- A:+)2 ^ ^ ^ 

1 



((P' - A;)-, - (P' - k)l){P'- - (P' - - fc„ 
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+ 



1 

iP'~ - iP' - k)-n - k~n){P- - P'- + {P' - k)-^ - (P - k)i 
{P' - fc)- ^ (P' - k) 



The third term in Eq. 



R 



, Fy^)^ is 



(53) 



27r 



d^k±dk+dk-e{k+ - p+)e{P'+ - k- 

(A;+)2(P+ - fc+)(P'+ - k+y 

1 



-Tr[0^{k-)] 



X 



X 



(fc- - fc- + ie){P- -k--iP- fc)- - ze)(P'- -k'- (P' - /c)- + le) 
1 



(A;- - A;^ + 2e)(P- - k- - {P - k) 



R 



It 



(54) 



In Eq. (^) the only pole of the integrand in the upper complex semi-plane of k is kf^^-^, 
given by Eq. (|I6]). By computing the residue, the result is 



-N 



d'^k^dk+0{k+ - P+)e{P'+ - V 

{k+y{p+ - k+){p'+ - k+y 
1 



'-Trp^P'- - (P' - k)-J] 



X 



(P'- - (P' - k)^ - k-J{P- - P'- + (P' - k)-„, - (P - k)-„ 

1 

(P'- - (P' - A:)-„ - k^){P- - P'- + (P' - A:)-„ - (P - Z.)^) 
The last term is Eq. (^91) , P^^^*^, is given by 



(55) 



-AT 



d'^k^dk+dk-e{k+ - p+)e{P'+ - k+) 
2ti^' J k+{p+ - k+y{p'+ - k+y 

1 



Tr[0+{k-)] 



X 



ik~ - fc- + ze)(P- -k--{P- k)-^ - it){P'- -kr- (P' - /c)- + le) 

1 

(P- -k--{P- k)-^ - It) iP'- -k-- (P' - k)-^ + le) • 



(56) 



In Eq. ([56| ) the poles /^^.g-j, Eq. (|1^), and /c^-^^, Eq. ([17|), contribute. By calculating the 
residues of the integrand of Eq. (|56|), we obtain 



-N 



d^k^dk+e{k+ - P+)e{P'+ - A;^ 

k+{p+ - k+y{p'+ - k+y 
1 



-Tr[O^P'- - {P' - k)-J] 



X 



X 



+ 



(P'- - (P' - fc)o-. - k-J (P- - P'- + (P' - A:)-„ - (P - k)-J 

1 



(P- - P'- + (P' - /.)-„ - (P - k)],) ((P' - A;)-„ - (P' - fc) 
- k)-^ ^ {P' - k)] 



R) 



(57) 
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FIGURE CAPTIONS 

Fig. 1. Light-front time-ordered diagrams for the current: (a) (Eq. (|^)) and (b) 
F^^) (Eq. (I)). 

Fig. 2. The pion form factor vs — g^. The contributions to F-^(g^), evaluated in the 
Breit frame with = \J —q^^ i.e. a = 90" (Eq. (P), are also shown. Sohd line: full result; 
dashed line: F^^\q^,a) (full result without the pair term, Eq. (|^)); dotted line: F^^^\q'^,a) 
(pair term, Eq. (||)); long-dashed line: Fy]„^^(g^,a) ( instantaneous-term contribution to 
Fy)(g^,a)); short-dashed line: F^^ 2ist{'f i ( instantaneous-term contribution to the pair 
term). Experimental data: Ref. |^ (full squares) , Ref. (full triangles) , Ref. |^ (empty 
squares), Ref. (empty circles) and Ref. (full circles). 

Fig. 3. Contributions to the pion form factor vs a (see Eq. (H)) for —q^ = 1 (GeV/c)^. 
Solid line: full result; dashed line: Fy\q'^,a) (full result without the pair term, Eq. (0)); 
dotted line: F'^^^\q^,a) (pair term, Eq. (|)); long-dashed line: Fy]„^j(g^, a) ( instantaneous- 
term contribution to F^^\q^^a))] short-dashed line: Fy£^j(g^,a) (instantaneous-term con- 
tribution to the pair term). 

Fig. 4. Comparison between the pion form factor calculations for (Eq. (^) 
and for F^^^) (Eq. (^), at a = 0°, 45° and 90°. Solid line: F^); dotted line: F^^^\ 
normalized to 1 at = 0. 

Fig. 5. Form factor for a weakly relativistic spin-0 system vs —q^. The contri- 
butions to the form factor, evaluated in the Breit frame with g"*" = a/— g^, i.e. a = 90° 
(Eq. (^)), are also shown. Solid line: full result; dashed line: Fy\q'^,a) (full result with- 
out the pair term, Eq. (^); dotted line: F^^^q'^, a) (pair term, Eq. (|D); long-dashed 
line: |Fy]„^^(g^, a)| ( instantaneous-term contribution to Fy)(g^,a)); short-dashed line: 
Fy£^((g^,a) (instantaneous-term contribution to the pair term). 

Fig. 6. Dimensionless product of the transverse momentum probability density (Eq. 
(P^)) times k1_ vs the ratio k±/m. Results for the present pion model (solid curve with dots), 
the pion model of Ref. [15] (dashed curve), a weakly relativistic system (solid curve) and the 
deuteron from the Reid soft-core potential (short-dashed curve). 

Fig. 7. Ratio between the pair term, calculated for a = 90°, and the full form factor 
vs — g^. Upper curves correspond to the pion. Solid line: covariant pion model, R^'^'^^^ (Eq. 
([33|); dashed line with dots: FFHD ratio, R^^^^^^ (Eq. (|35D); short-dashed line: kinematical 
estimate,i?(^^^) (Eq. (g)). Lower curves correspond to a weakly relativistic system at the 
deuteron scale. Solid line: covariant result; dashed line: kinematical estimate. 
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